607 
| VOLUME 81 SEPARATE No. 607 


PROCEEDINGS 


AMERICAN SOCIETY 
OF 


CiVIL ENGINEERS 
JANUARY, 1955 


LATERAL BUCKLING OF ECCENTRICALLY 
LOADED I-COLUMNS 


by Mario G. Salvadori, M. ASCE 


ENGINEERING MECHANICS 
DIVISION 
{Discussion open until May 1, 1955} 


Copyright 1955 by the American Society or Civit ENGINEERS 
Printed in the United States of America 


Headquarters of the Society 
33 W. 39th Se. 
New York 18, N. Y. 


PRICE $0.50 PER COPY 


nT 


THIS PAPER 


--represents an effort by the Society to deliver 
technical data direct from the author to the 
reader with the greatest possible speed. Tothis 
end, it has had none of the usual editing required 
in more formal publication procedures. 


Readers are invited to submit discussion apply - 
ing to current papers. For this paper the final 
date on which a discussion should reach the 
Manager of Technical Publications appears on 
the front cover. 


Those who are planning papers or discussions 
for “Proceedings” will expedite Division and 
Committee action measurably by first studying 


“Publication Procedure for Technical Papers” 
(Proceedings — Separate No. 290). For free 
copies of this Separate—describing style, con- 
tent, and format—address the Manager, Techni- 
cal Publications, ASCE. 


Reprints from this publication may be made on 
condition that the full title of paper, name of 
author, page reference, and date of publication 
by the Society are given. 


The Society isnot responsible for any statement 
made or opinion expressed in its publications. 


This paper was published at 1745 S. State Street, 
Ann Arbor, Mich., by the American Society of 
Civil Engineers. Editorial and General Offices 
are at 33 West Thirty-ninth Street, New York 18, 
N. Y. 


| 


LATERAL BUCKLING OF ECCENTRICALLY LOADED I-COLUMNS 


Mario G. Salvadori,! M. ASCE 


SYNOPSIS 


The critical combination of thrust and unequal end- moments or eccentrici- 
ties (interaction curves) for the lateral buckling of columns elastically re- 
strained against rotation about the principal axes at the end sections, and 
prevented from twisting about the column axis at the end sections are derived 
for various values of the Timoshenko parameter L2/a2 and of the end- 
moments ratio. 

It is shown that for most practical cases these critical combinations are 
easily obtainable as soon as the critical value of equal end-moments without 
thrust is known for columns free to rotate at their ends (simply supported 
columns). 


INTRODUCTION 


The solution of the lateral stability problem for I-columns under eccentric 
thrust, i.e., under a combination of axial thrust and end- moments, consists in 
the determination of so-called interaction curves from which the critical 
value of the thrust may be derived for given end- moments or the critical 
value of the end-moments for a given thrust. 

The shape of the interaction curves depends essentially on: 1) the 
geometric and strength characteristics of the I-column and, in particular, 
on the value of the Timoshenko parameter L2/a2 ; 2) the ratio r of the end- 
moments, or eccentricities; 3) the boundary conditions at the column ends. 

Interaction curves for I-columns with equal end- moments were obtained 
by Johnston,(1)2 for columns of rectangular cross-section with unequal end- 
moments by DiMaggio,(2) for I-columns with unequal end- moments by 
Salvadori.(3) All these solutions considered “simply supported” columns, 
whose ends are free to rotate around their principal axes at the ends and are 
prevented from twisting around the column axis there. 

In the present paper interaction curves are obtained for I-columns whose 
ends are free to rotate in the plane of the web (strong plane), but are either 
elastically restrained or totally fixed in the plane of the flanges (weak plane) 
at the ends, while they are prevented from twisting around the column axis 
there. The consideration of unequal end- moments in the plane of the web, 
which may be due to continuity, allows the use of these interaction curves for 
columns elastically restrained in both the weak and the strong planes. 

The treatment is limited to stresses within the elastic limit. 


1. Prof. of Civ. Eng., Columbia Univ., New York, N. Y. 
2. Numbers in parentheses refer to the bibliography at the end of the paper. 
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List of Symbols 


on Wn 


column length 

aL = reduced column length 

column depth 

minimum flexural rigidity of the column 

torsional rigidity of the column without thrust 

flexural rigidity of one flange in its own plane 

bending moment at top of column in the plane of the web 
bending moment at top of column in the plane of the web without 
thrust 

bending moment at the bottom of column in the plane of the web 
M,/M, = end-moment ratio 


MoL/ VB C = bending moment coefficient 

bending moment coefficient without thrust 

bending moment coefficient for equal end-moments without thrust 

My {/ VBC = aK, = bending moment coefficient in terms of 
reduced length 

a Ky 1 = reduced length bending moment coefficient for equal 

end- moments 

Ky, r/Ko 1 = magnification coefficient for unequal end- moments 

minimum, safe value of 

(2L2c)/ (h2pD) = Timoshenko parameter for I-beam lateral buckling 

a*L2/a2 = Timoshenko parameter in terms of reduced length 

axial thrust 

nm 2p/ (2 = minimum Euler thrust 

P/P., = thrust coefficient 

critical value of equal end-moments without thrust 

Mp ;/P = top and bottom eccentricities 

Kr/Ko,r = Mp/Mg 2 = ratio of top moments with and without thrust 


eo VB/CL? = top eccentricity coefficient 
CA/Ip = critical torsional thrust 


= cross-sectional area of column 


polar moment of inertia of column cross-section 
C(1 - P/Py) = reduced torsional rigidity 
P.,/ Py = ratio of minimum Euler thrust to torsional critical 
thrust 
Euler thrust in the plane of the web 
total energy of buckled column 
lateral displacement of buckled column in the plane of minimum 
rigidity 
angular rotation of buckled column 
abscissa along column axis with origin at bottom of column 
My [r+ (1-r)z/L] = bending moment in the plane of the column 
web 
8’, 8’; ... = derivatives of u and 8 with respect to z 
(u/h) /2B/D = non-dimensional lateral displacement 
functions of z for the expansion of 7 and f 
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Mj = 
Tr = 
K, = 
Kor 
Koi 
ky = 
Ko 1 
= 
= 
L2/a2_ = 
(2/a2 a 
Pp = 
Por 
p s 
Mcr 
€2,1 
Kr 
Pr = 
A 
Ip = 
H s 
Pe = 
Vv 
u = 
B 
vA = 
M(z) = 
= 


Lateral Buckling under Equal End- Moments 


A “simply supported” column of length L, depth h, minimum flexural 
rigidity B (in the plane of the flanges), torsional rigidity C, and flange 
rigidity D buckles paren ' under the action of two equal end- moments M, 
in the strong plane when: (4) 


n_/BC \° 
i, 1 + (1) 


Introducing the moment coefficient:3 


Koa =ML/ , (3) 


its critical value is given by: 
n \2 
Koa = yi + (4) 


When the I-column is completely fixed in the weak plane at the ends, the 
critical value of the end-moments is given by: 


(5) 


and hence: 


. (6) 


K ] Qn 1+ 


The critical value of K, ; may be expressed by means of a single formula 
if we introduce the concept of reduced length f, i.e., of the length between 


inflection points when the same column buckles under thrust (Euler buckling). 
Letting 


L=eat (1/2 € (7) 


the reduction coefficient @ equals 1 for columns simply supported in the weak 
plane, 1/2 for columns fixed in the weak plane, and has values between 1/2 
and 1 for columns elastically restrained in the weak plane. By means of { 
the expression for Ko 1 becomes, whatever the conditions of restraint in the 


weak plane: 


3. The first subscript, 0, indicates the absence of thrust; the second, 1, the 
equality of end-moments, whose ratio r equals one. 
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where: 
a“ D 
— 


= MA / BC = 1+ 


and the critical value o. the moments is given by: 


cr O,l L 


An analogous result wis obtained by Mr. G. Winter, M. ASCE, for beams of 
rectangular cross-sec'ion.(7) Figure 1 gives the values of Ky versus 
2/a , which also appvar in the row r = 1 of Table I, where { = L and 


hence ky = K, 


Lateral Buckling under Unequal End- Moments 
Indicating by: 


r=M (10) 


0,1/Mo,2 


the ratio of the moment M,, , at the bottom of the column to the moment M, 2 
at the top of the column, the critical values of M, 2 and of Mp j= rMp 2 can 


be obtained by the Rayleigh-Ritz method as shown in the Appendix (Sec. 10). 
Tables I and II give: 


= My (11) 


for columns simply supported in the weak plane (taken from ref.(3)) and for 
columns fixéd in the weak plane, respectively, in terms of L2/a2 and r. 

By means of the coefficients K, , of Tables I and Il it is easy to compute 
in terms of L2/a? the ratios: : 


of the critical moment Mo, 2 at the top of the column to the critical value 

Mo 2 ] r = 1 of the same moment when Mo 2 = Mo,1(r = 1) for the two ex- 
treme conditions of simple support and fixity in the weak plane. Table III 
shows that for all practical purposes this ratio is independent of L2/a2 for 
L2/a2 £ 40 and has the same value for both extreme conditions of support 
for r= 1/2 and r=0. For r = -1/2 andr = -1 this ratio is slightly sensitive 
to the values of L2/a2 and to the conditions of ee Hence, it may be 
conjectured that for all values of r and for L2/a2 © 40 the ratio H, will be 


comprised between the extreme values corresponding to the conditions of 
simple support and of fixity (Fig. 2). 

Noticing that in most practical cases L2/a2 ¢ 40 and that errors of up to 
10% are easily permissible if on the safe side, a simplified approximate 
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evaluation of the critical moments may be obtained by means of the minimum, 


safe, values H, of HM, given in Table IV or Fig. 2. 
By means of Eqs. (8) and (11) the critical value of My 2 may be obtained, 
whatever the degree of restraint in the weak plane, by means of the formula: 


yBC 


aL 


where k, 1 is given by Fig. 1. It is seen by comparison of Eq. (13) with 
Eq. (11) that: 


(14) 


As clearly shown by Eq. (14), the introduction of the (safe) ratios 4 r and 
of the reduced length f, reduces the solution of the lateral buckling problem 
— unequal end-moments, whatever the end restraints and the value of 
L2/ a2, to the known solution of the buckling problem of a simply supported 
column in pure bending. 

For example, an elastically restrained column of reduced length 
{ = 0.7 L (a@ = 0.7) with a moment Mp 2 at the top and no moment at the 


——- (Mo, 1= TM, 2 = 0, i.e., r = 0), for which L2/a2 = 32, has an 


(?/a2 = (0.7 L)*/a2 = 16. Hence, from Fig. 1 or Table I, ky ; = 4.00; from 
Fig. 2 or Table IV HZ, = 1.82, and by Eq. (13): 


_ 
Moo = 1.82 x 4.00 = 10. 4 /BC /L. (a) 


A more accurate evaluation of My 2 may always be obtained by means of the 
values 4, of Table III, but is seldom warranted in practice. 


Lateral Buckling under Equal End-Moments and Thrust 


Indicating by: 


2 
(15) 


cr (x L)* 


the minimum critical value of the thrust in the absence of end-moments, and 


by: 
f3C 


aL 
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(16) 

| 


the critical value of equal end-moments in the absence of thrust, the critical 
combination of thrust P and equal end-moments M, that is the expression for 
the interaction curve P versus M, is given by: 


2 
+ (M/M,_.) = (17) 


Letting: 


k, = M4/ 


we obtain from Eq. (17): 


For p = 0, i.e., no thrust, kK; = Ko 1 and the ratio: 
/k = fl-p (21) 
1° O,1 


is independent of L2/a2 and @. The curve ky/Ko 1 versus p is identical with 
the curve K,/K, 1 for a = 1 and appears labeled r = 1 in Fig. 3 (thin curve). 


By means of the values Ko 1 = Ko1 of Fig. 1 the critical value of equal 
end- moments with thrust is given by: 


(22) 


When the moments M are due to the eccentricity e of the thrust P, the 
critical value of e may be computed from the equation: 


e = M/P (23) 


where M is given by Eq. (22). 
Vice versa, Eq. (17) may be solved for P in terms of e, giving: 


2 


2 2 
PL 
n 
and: 
0,1 aL 
2e 
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Lateral Buckling under Unequal End- Moments and Thrust 


The critical combinations of thrust P and of top moment My for various 
values of the ratio r = M;/Mo and of L2/ a2 may be computed by the methods 


explained in the Appendix for conditions of simple support {computed from 
ref.(3), and for conditions of fixity in terms of the non-dimensional ratios 


p = P/Per (Eq. (18)) and: 


(25) 


= ypc . 


Remembering the definition of Ko r (Eq. (11)), it is also seen that: 


Table V gives the largest percentage difference in absolute value between 
the ratios Kk. for simply supported and fixed ends, i.e., the values: 


- & Je 
Ok,, = max. x 100. 


It is seen that A K, is less than 3% for all values of L2/a2 and of r, except for 
r=-l. Hence, for practical purposes, the interaction curves Ky versus p 


are independent of L?/a2 and of the conditions of restraint in the weak plane, 
except for r = -1. 


Figure 3 gives the graphs (thin curves) of the average value of K, versus 
pforr=1,r=1/2,r=0, andr=-1/2. These graphs can be used whatever 
the value of L2/a2 and the end restraints. For r = -1, Fig. 3 gives 3 separate 
graphs for ™ a supported ends and 2 separate graphs for fixed ends in 


terms of L2/a4. 

The discontinuity in the slope of the curves r = -1 of Fig. 3 corresponds to 
the meeting point of two separate branches of the interaction curves. The 
branch valid for lower values of x, represents a buckling deformation in 


which the flanges are displaced in the same direction; the branch valid for 
x, greater than 0.8 corresponds to flange displacements in opposite direc- 


tions. This last branch if prolonged to meet the p axis, would cross it ata 
value corresponding to the second Euler load, i.e., at p = 4. 
By means of the ratios K, of Eq. (27) the top moment Mz may be written: 


(28) 


or by means of Eq. (13): 


To illustrate the use of Fig. 3 consider a column restrained in the weak 
plane so as to have a reduced length ( = 0.7 L(a@ = 0.7), with L2/a2 32, and 
loaded by a thrust P = 0.6 P... Assume the restraint conditions in the strong 
plane to be such that r = 0. With p = 0.6 andr = 0 we obtain from Fig. 3 
Ky = 0.66. For r = 0 Table IV (or Fig. 2) gives H, = 1.82, for 
(2/a2 = (0.7)2 x 32 = 16 Table I (or Fig. 1) gives ky 4 = 4-00. Hence, by 
Eq. (30): 


3 0.66 x 1.82 x 4.00 /BC BC 


The corresponding eccentricity eg» of the thrust at the top of the column is 
given by: 


=M./P= yBC/L cx? 
e, M,/P M,/(0.6 0.57 / CL 


0.6n°B/(0.7L)* 


To solve the instability problem for a given eccentricity: 


€, CL°/B , 


notice that M, = P ey = p P,, eg and that by Eq. (26) Mz = Ky VBC/L. 
Hence: 


L)* 2 B 


from which: 


or, dividing both sides of this equation by Ky, ;: 


BC 
2) 
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(33) 


The intersection of the straight line p versus x, of Eq. (33) with the p versus 
Ky graph of Fig. 3 gives the critical value of p and hence of P. 
For example, for eg = 0.4 VCL2/B, i.e., €9 = 0.4, @ = 0.7, L?/a? = 32, 


r = 0, we have (2/a2 = 16 and by Eq. (14), Table IV, and Fig. 1 K, , = 


ko 1 = 1.82 x 4.00/0.7 = 14.8. Hence: 


0.7 2 x 14.8 


"9.87 x 0. = 1.83. 


The line through the origin of slope equal to 1.83 intersects the dashed curve 
r = 0 of Fig. 3 at p = 0.84; therefore, the critical value of the thrust equals: 


P = 0.84 «°B/(0.7L) = 1.20 (a) 


Reduction of Torsional Rigidity due to Thrust# 


A column buckles in pure torsion under a thrust P,, uniformly distributed 
over its cross-section when:(5) 


CA/ (34) 


where A is the area of the cross-section and i) its polar moment of inertia. 


The equivalent or reduced torsional rigidity H of the column under a thrust 
P is given by:(5,6) 


H=C(1- P/P,) =C(l-p (35) 


In order to take into account the reduction of torsional rigidity in the evaluation 
of the interaction curves, the moment coefficient K,. is defined in terms of H 
rather than in terms of C: 


= /BH = (M,L/ /BC) /C/H = K, . (36) 


4. This section is based on a suggestion of Mr. J. W. Clark. 


aK, 
Es 


4 


Noticing that for P = 0 the reduced rigidity H, equals C and indicating by 
My 2 the value of M2 in the absence of thrust: 


= z = 
M, ob//BH, = M, = K, (37) 


Hence the ratio of moments M,/Mo 2 becomes: 


The ratio: 


I 


Y = = (39) 


T cA 


varies, practically, between zero and one. When the reduction of torsional 
rigidity is negligible, y = 0 and the thin curves of Fig. 3 remain unchanged. 
When the torsional rigidity is greatly reduced by the thrust, y= 1 and the 
ordinates kK. of these curves must be multiplied by ¥1- p. The thick curves 
of Fig. 3 give ky versus p for y= 1 and are on the safe side for all practical 
cases. For intermediate values of y new curves may be easily computed. 

To illustrate the influence of the reduction of torsional rigidity on the value 
of the critical moments, consider the column of the example of Sec. 5. As- 
suming the highest possible value of y(y = 1) and p = 0.6, the moment My 
of Eq. (b) becomes: 


M, = yl - 0.6 x 6.87 /BC /L= 4.36 yBC /L. 


Influence of Deflections in the Plane of the Web 


The interaction curves of the preceeding Sections are obtained under the 
assumption that the flexural rigidity of the column in the plane of the web is 
much greater than the rigidity in the plane of the flanges, so that bending 
deflections in the plane of the web may be neglected. 

For the case of equal end-moments the influence of the deflections in the 
plane of the web may be taken into account by multiplying the critical bending 
moments M by the factor 1 - P/Pg, where Pg is the Euler thrust in the plane 


of the web. Thus the curves r = 1 of Fig. 3 may be used provided: 


K - (40) 


be substituted for Kr: 


(38) 
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The influence of the deflections in the plane of the web decreases as r 
approaches -1. Hence this influence may be taken empirically into account, 
with sufficient accuracy, by substituting for x, the ordinates: 


- P (r + 1) (41) 


In the examples of Secs. 5 and 6, assuming a ratio Per/PE = Imin/Imax 
= 0.2, with r = 0 and p = 0.6 the critical moment becomes: 


M, = 4.36 (1-0.5 x 0.6 x0.2) /L = 4.10 /L 


CONC LUSIONS 


The results of the preceding Sections prove that for all practical purposes 
the critical combinations of thrust and end- moments of I-columns elastically 
restrained in bending about the principal axes of their end-sections and fixed 
against twisting there, may be easily obtained once the critical value of equal 
end- moments for a simply supported column is known. By means of the 
graphs and tables presented in this paper the designer may evaluate critical 
thrusts and moments (or eccentricities) with a minimum of computations. 

The most significant result achieved by this numerical research is perhaps 
the discovery of the practically linear relation between p and x, valid for 
most practical cases when the reduction of torsional rigidity is taken into 
account. 

Although it must be remembered that an extension of the present results 
to the inelastic range is essential to all engineering applications, it is hoped 
that their simplicity will encourage specification writers to the adoption of 
more accurate column formulas in order to achieve greater economy and 
safety in steel structures. 


Mathematical Appendix 


The values of the moment coefficients: 


K. = (42) 


and of the ratios x, = K,/K, r Were obtained by the Rayleigh-Ritz method 


starting from the expression for the total energy V of an I-column buckled 
laterally under the action of thrusts P and end- moments: 3 


2 


+ + [up u"dz - (43) 


(see “List of Symbols” for meaning u, 8, z and M). 
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The boundary conditions for columns simply supported in the weak plane are: 


u(O) = u"(0O) = u(L) = u"(L) = 


(44) 


(0) = p"(0) =A(L) = p'(L) = 


For columns fixed at the ends in the weak plane the boundary conditions 
are: 


u(O) = u'(0) = u(L) = u'(L) = 0 


(45) 


p(o) = p'(o) = = = 


The conditions of Eqs. (44) are satisfied by taking: 


sin nx z/L ; 6 (2) = 2 


u(z) 1 ne z/L . (46) 


1 
The corresponding Rayleigh-Ritz calculations are summarized in ref.(3) 
To satisfy the conditions of Eqs. (45), let: 


D 
u(z) =h (z) (47) 
and expand 7 and # in the following series: 


(2) =, +2, © n (48) 


where: 
(z) = (1 - cos 2nxa z/L) (50) 
y (2) = [2z/L -l1-sin7 (22/L - 1)/sin (51) 
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and vp are roots of the transcendental equation: 


tan Y = (52) 


The choice of these functions greatly simplifies the calculations, since 
they are orthogonal together with their first and second derivatives in almost 
all the combinations appearing in the integrals to be evaluated in minimizing 
the total energy V. 

In order to obtain an accuracy of 1% or better the cases r = 1/2, r = 0, 
r=-1/2,and r=-1 were solved taking into account one term of each series 
and from the corresponding values K, were computed the ratios x, = K,/Ko,r- 
The corresponding determinantal equations are of the 4th order. The case 
r = -1 was then solved with two terms of each series, leading to determinantal 
equations of the 8th order. The new K, differed substantially from those of 


the first approximation, but the new ratios K, were instead practically 


identical with those obtained by means of the first approximation. Hence the 
ratios x, of the first approximation were considered correct within the re- 
quired accuracy and only the Ky , were computed again by means of 2 terms 
of each series for r = 1/2, r = 0, andr = - 1/2, and with 2 and 3 terms of each 
series for r = -1. In this last case the determinantal equations of order 8 
and 12 split, for p = 0, into two separate equations of order 4 and 6, whose 
roots correspond to the two branches of the curves r = -1 of Fig. 3. The 
intersection of the two branches was also determined in order to draw the 
graphs with greater accuracy. 

The extremely complicated calculations were carried out by the staff of the 
Istituto Nazionale per le Applicazioni del Calcolo, the mathematical laboratory 
of the Italian National Research Council directed by Professor Mauro Picone. 
Professor Aldo Ghizzetti, Vice-Director of the Institute, and Dr. W. Gross 
supervised the calculations. Mrs. E. D. Broadbent, B.Sc., Dr. C. Mengotti 
Marzolla, and Mr. A. Perugini were in charge of the calculations, Mr. V. 
Ceccarini participated in the calculations and drew the graphs. 

Iam greatly indebted to the Istituto Nazionale per le Applicazioni del 
Calcolo for this accurate and laborious numerical research, which was 
partially supported by a grant of the Column Research Council of the 
Engineering Foundation, under the chairmanship of Mr. S. Hardesty, and 
which was sponsored by its Committee E, headed by Dr. H. V. Hill. 
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